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ABSTRACT 



Aims. We constrain a hypothetical variation in the fundamental physical constants over the course of cosmic time. 
Methods. We use unique observations of the CO(7-6) rotational line and the [Cl] 3 P2 — 3 Pi fine-structure line towards 
a lensed galaxy at redshift z = 5.2 to constrain temporal variations in the constant F — a 2 / fi, where fj, is the electron- 
to-proton mass ratio and a is the fine-structure constant. The relative change in F between 2 = and z = 5.2, 
AF/F= (F b s — fi a b)/Fi a b, is estimated from the radial velocity offset, AV = V I0 t — V&, between the rotational 
transitions in carbon monoxide and the fine-structure transition in atomic carbon. 
Results. We find a conservative value AV = (1 ± 5) km s _1 (la C.L 
gives AF/F < 2 x 10 
(look-back time t z o.7 - 

limit on Afi/fi < 2 x 10~ 7 at z = 5.2 (t Z 5.2 — 12.9 Gyr). From both constraints on AF/F and Afi/fj, , one obtains 
for the relative change in a the estimate Aa/a < 8 x 10 -6 , which is at present the tightest limit on Aa/a at early 
cosmological epochs. 



which when interpreted in terms of AF/F 
5 . Independent methods restrict the /i- variations at the level of Afj,/fi < 1 x 10~ 7 at z — 0.7 
6.4 Gyr). Assuming that temporal variations in fi, if any, are linear, this leads to an upper 
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Techniques: radial velocities — Elementary particles 
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1. Introduction 

Space-time variations n the dimensionless physical con- 
stants such as the fine-structure constant, a = e 2 /(hc), 
and the electron-to-proton mass ratio, fj, — m c /m p , are 
predicted within the framework of grand unification the- 
ories, multidimensional theories, and scalar field theories 
(for a review, see Uzan 2011). These constants mediate the 
strengths of fundamental forces. In essence, a is the cou- 
pling constant of the electromagnetic interaction and m G is 
related to the vacuum expectation value of the Higgs field, 
i.e., to the scale of the weak nuclear force, whereas m p is 
proportional to the quantum chromodynamic scale, Aqcd ■ 
Hence, /i probes the ratio of the weak to strong forces of 
nature. If detected, a variation in any of these constants 
would be a manifestation of the need to develop new theo- 
ries since the Standard Model of particle physics does not 
predict any variations in fundamental constants. 

Experimentally, the variations in the dimensionless con- 
stants can be measured by analyzing the relative positions 
of atomic and molecular transitions measured in both lab- 
oratory and astronomical objects. At present, the most ac- 
curate laboratory results on temporal a- and /^-variations 
are a/a = (1.6 ±2.3) x lO" 17 yr" 1 (Rosenband et al. 2008) 
and jjl/jjl = (1.6 ± 1.7) x 10~ 15 yr" 1 (Blatt et al. 2008). 
For values of a(t) and n(t) that are linearly dependent on 
cosmic time t, these limits correspond to constraints at red- 



shift z ~ 2 (look-back time t z2 ~ 10 Gyr) of |Aa/a| < 0.4 
ppm (lppm = 10~ 6 ) and |A/t//t| < 30 ppm, where Aa/a 
(or A/t//t) is a fractional change in a between a reference 
value cti and a value a 2 measured in a cosmic object given 
by Aa/a — (a 2 — a\)/a\. 

The most accurate astronomical measurements of the 
cosmological [i- variation are performed in the radio range, 
which constrain Au./u. at the level of 0.3 ppm at z = 0.89 
(Ellingsen et al. 2012) and 0.1 ppm at z = 0.69 (Kanekar 
2011). The lowest spectroscopic limits on a are \Aa/a\ < 3 
ppm at z = 0.77 obtained from radio observations (Kanekar 
et al. 2012) and the same order of magnitude |Aa/a| < 3 
ppm at z ~ 2 — 3 derived from optical spectra of quasars 
(Agafonova et al. 2011; Molaro et al. 2008; Srianand et 
al. 2007; Quast et al. 2004). To complement these results, 
the measurements of Murphy et al. (2003) claim detection 
of an a-variation in about 140 absorption systems from 
quasar spectra observed with the Keck telescope Aa/a 
= —5.4 ±1.2 ppm in the range 0.2 < z < 3.7. However, 
Griest et al. (2010) showed that the wavelength calibration 
of these Keck spectra was affected by systematic and un- 
explained errors of ~ 500 m s _1 , which transformed into 
an error of o~& a / a ~ 30 ppm. The Keck result was not 
confirmed by the analysis of the VLT spectra of quasars ac- 
quired from the southern hemisphere by Webb et al. (2011), 
who found a positive signal of Aa/a = 6.1 ± 2.0 ppm at 
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z > 1.8 in about 150 systems, and proposed a model of a 
cosmic dipole to reconcile their Keck and VLT findings. A 
tentative detection of the change in a of Aa/a = — 3.1±1.2 
ppm was also reported by Kanekar et al. (2010) based on 
observations of OH 18 cm lines at z = 0.25. 

At earlier cosmological epochs z > 4, only a handful 
of spectral measurements have been able to provide direct 
constraints on the temporal variations in a and \i. Some 
time ago, we suggested a few methods to probe these vari- 
ations at high redshifts by using far-infrared, sub-mm, and 
mm transitions in atoms, ions, and molecules usually ob- 
served in Galactic and extragalactic sources (Kozlov et al. 
2008; Levshakov et al. 2008; Levshakov et al. 2010a). The 
advantage of fine-structure (FS) transitions is that they 
have the sensitivity to an a- variation that is about 30 times 
higher than that of UV transitions employed in optical spec- 
troscopy (Levshakov et al. 2008, hereafter L08): if Aa/a 
7^ 0, then the frequency shift Auj/lu{ s = (u) z ,ts — <^fs)/wfs = 
2Aa/a. To take advantage of the high sensitivity of the 
FS transitions in differential measurements of a, we sug- 
gested to take as a reference the low lying rotational lines 
of carbon monoxide CO since the rotational frequencies of 
light molecules are independent of a but sensitive to the 
value of /i, Au;/u! ro t — Afi/fi. In case of variations in a 
and/or \i over cosmic time, one should observe the ap- 
parent redshifts of the FS and rotational lines differ of 
Az/(1 + z) = AV/c = AF/F, where F = a 2 //i, AV 
is the difference of the apparent radial velocities, and c 
is the speed of light (L08). The implementation of this 
method has resulted in AF/F < 100 ppm at z = 6.42, 
and AF/F < 150 ppm at z = 4.69 (L08), and AF/F < 85 
ppm over the redshift interval z = 2.3 — 4.1 (Curran et al. 
2011). 

In the present Letter, we show that the [Ci]/CO emis- 
sion discovered by Combes et al. (2012, hereafter CRR) to- 
wards the lensed galaxy HLSJ091828.6+514223 at z = 5.2 
allows us to significantly improve this limit. We note that 
the look-back time is t z $.2 = 12.9 Gyr for the cosmological 
parameters H : ft b : il c : il A = 70.5 : 0.046 : 0.228 : 0.726 
(Hinshaw et al. 2009), which is about 94% of the total age 
of the Universe. 



2. Observations and results 

Details on observations of the z — 5.2429 lensed galaxy 
HLSJ091828.6+514223 are given in CRR. Here we con- 
centrate on the CO(7-6)/[Ci](2-l) data since both lines 
were observed simultaneously with the same 2mm receiver 
at the IRAM-interferometer (PdBI). The spectra were re- 
duced with the channel width of A c h = 46.4 km s _1 , and 
the resulting noise rms was approximately 0.95 mJy per 
channel. 

The observed spectra of CO(J = 7-6) and [C i](J = 2-1) 
are shown in Fig. [1] We clearly detect two subcomponents 
separated by ~ 510 km s _1 and provide in Table 2 in CRR 
the parameters for a wider and stronger 'red' component of 
linewidth (FWHM) w r ~ 540 km s _1 , and a narrow and 
weaker 'blue' component of Wb — 150 km s _1 . We note that 
the CO(7-6) and [Ci](2-1) lines, thanks to their close rest- 
frame frequencies (806651.806 MHz and 809341.97 MHz, 
respectively), were observed within one bandwidth, thus 
their relative positions are free from possible instrumental 
systematics. 



In CRR, the profiles of [Ci](2-1) and CO (7-6) were fit- 
ted by two simple Gaussians (for the blue and red compo- 
nents) centered at V& = — 530±9 km s _1 and 20±14 km s _1 
for the former line and at V r = — 510±30 km s _1 and 4±10 
km s _1 for the latter (Table 2 in CRR), i.e., the velocity off- 
set between them is AV = V TOt — Vf s = 20±31 km s _1 (blue) 
and — 16 ± 17 km s _1 (red). We note that the precision of 
a single Gaussian line center is a — 0.69 • rms ■ i/Ach . w/Sq 
(Landman et al. 1982), where 5*o is the peak intensity and 
w the linewidth. For the parameters w ~ 150 km s _1 , 
A c h ~ 46 km s~ x rms ~ 1 mJy, and So ~ 10 mJy (cf. 
Fig. [J), one finds a ~ 6 km s _1 , which is slightly smaller 
than the error oy h quoted above. 

In the present study, we employ a different - model- 
free - method to measure AV. We adopt our approach be- 
cause despite their complexity, the CO(7-6) and [Ci](2-1) 
lines have almost identical shapes (the scaling factor be- 
tween their intensities is ~ 1.6, as indicated in Table 2 in 
CRR). The procedure ("sliding distance") is realized in the 
following way: one of the chosen profiles (e.g., T\ for CO) 
is fixed and the other (T 2 for [C 1] ) is shifted relative to T\ 
sequentially in small steps, 5v, within the velocity interval 
[V, V2]. At each step, a value of O 2 - the sum of squares of 
the intensity differences over n points within the line pro- 
file - is to be calculated Q 2 (5v) = £"=ipi,j - T 2j (8v)} 2 , 

where O is a kind of distance between the curves, and O 2 
depends parabolically on Sv that has its minimum value 
when the profiles are aligned. The velocity at this mini- 
mum is taken as the offset, AV, between the two profiles. 

In the rest frame, the Doppler velocity shift between 
the frequencies of [Ci](2-1) and CO(7-6) is AV = -999.8 
km s _1 . From Fig. [TJ it can be seen that there is a weak 
component at V ~ -1800 km s _1 in the [Ci](2-1) profilcQ 
that may be present in the CO (7-6) profile as well, thus 
distorting part of the [Cl](2-1) profile around V ~ —800 
km s _1 . To eliminate the influence of this blending on the 
minimum of O 2 , the intensity differences were calculated 
only in the velocity window —650 < V < 280 km s _1 
(marked by the vertical lines in Fig. [1]). In addition, to 
ensure approximately the same intensities in both profiles, 
the intensity of [Ci](2-1) was multiplied by 1.6 (Fig. H]). 

The function O 2 takes a minimum value at AV = 0.8 
km s _1 (red line in Fig. [3]), i.e., the [C 1] (2-1) line is shifted 
by —0.8 km s" 1 relative to CO(7-6). The error in this offset 
can be estimated by a bootstrap method, where artificial 
data samples are created based on the statistical properties 
of the reference data set. In the present context, this means 
that we have to create spectral profiles that are statistically 
equivalent to the observed ones. In Fig.[TJ we find both blue- 
and redward of the emission profiles the extended parts of 
the noise spectrum. The statistical analysis of these parts of 
the spectrum shows that the noise fluctuations are purely 
Gaussian (i.e., uncorrelated, with zero skewness and excess 
kurtosis) with a dispersion of 0.95 mJy. In a subsequent 
step, we obtain an emission- line template, i.e., a noise- free 
emission-line spectrum. Since the fitting of the Gaussian 
profiles may be uncertain because of the above-mentioned 
blending, the template was prepared by filtering the original 
data. We used a symmetrical Savitzky-Golay filter, which 
does not introduce any shifts into the filtered data (Savitzky 



1 This negative velocity wing in [Cl](2-1) could be related to 
the gas outflows (galactic winds) often seen in ultra-luminous 
infrared galaxies (CRR). 
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Fig. 1. The CO(7-6) and [Cl](2-1) emission lines 
(blue color histogram) towards the lensed galaxy 
HLSJ091828.6+514223. The velocity scale is given relative 
to z = 5.2429. The red component of the CO(7-6) line 
is centered at zero velocity. Its blue component is seen at 
V = —510 km s _1 . The red and blue components of [C i](2- 
1) are shifted with respect to CO(7-6) at —999.8 km s _1 
owing to their different rest-frame frequencies. Two verti- 
cal lines indicate the window used in the bootstrap calcu- 
lations. The residuals between the original and smoothed 
data (red color histogram) are shown at the bottom. 



& Golay 1964). The filtered profile is shown by the red 
color histogram in Fig. [TJ The size of the filtering window 
(n = 7) was chosen from the condition that the dispersion 
in the residuals (differences between observed and filtered 
data) in the velocity range of the observed emission should 
be equal to the dispersion in the noise in the emission-free 
regions (Fig. [1]). Artificial emission spectra (1000-10000 re- 
alizations) were generated by adding the noise fluctuations 
taken from the normal distribution of zero mean and disper- 
sion of 0.95 mjy to the template and the minimization of O 2 
was performed as described above. Examples of the <d 2 (Sv) 
curves for single realizations are shown in Fig. [3] The re- 
sulting distribution of the velocity offsets AV between the 
CO (7-6) and [Ci](2-1) profiles is plotted in Fig. [4] It is 
purely normal with the mean of 0.8 km s _1 and dispersion 
of 4.6 km s _1 . 

In calculations, different windows (indicated by verti- 
cal lines in Fig. [JJ were used to estimate the influence of 
the possible blending in the pair [Cl]/CO on the value of 
the velocity shift. The dispersion of the AV distribution 
was found to be very stable, varying by only a few percent, 
whereas the mean of the distribution slightly shifted from 
the above value of 0.8 km s" 1 when [Vi, V 2 ] = [-650, 280] 
km s- 1 to 0.07 km s" 1 at [VuVi] = [-650,0] km s" 1 . 
Thus, as a conservative estimate the value of AV = 1 ± 5 
km s _1 (lcr C.L.) can be taken for the velocity offset be- 
tween the CO(7-6) and [Ci](2-1) lines. If one interprets this 
offset in terms of AF/F , then AF/F < 17 ppm. Taking 
into account that Afi/fi < 0.1 ppm at t z0 , 7 = 6.4 Gyr and 
that the linearly extrapolated limit at ^5.2 = 12.9 Gyr, 
An/ n < 0.2 ppm, is two orders of magnitude lower than the 
limit on AF/F , one finds a constraint on the a-variation 
at z = 5.2 of Aa/a < 8 ppm. 



-2000 -1500 -1000 -500 500 1000 

Relative Velocity, km s _1 

Fig. 2. Same as Fig. [TJ but the original spectrum is shifted 
to +999.8 km s _1 (red color histogram) and its intensity is 
zoomed by a factor 1.6 to be comparable with the CO(7-6) 
spectrum within the window marked by two vertical lines. 
This window is used in the bootstrap calculations. 



3. Discussion and conclusions 

By minimizing O 2 , we have estimated the velocity offset 
AV between the CO(7-6) and [Ci](2-1) profiles with the 
precision of 5 km s _1 , which has provided an upper limit 
to the F variability at z = 5.2 that is an order of mag- 
nitude lower then the limit found in a previous survey of 
the redshifted pairs of [Ci]/CO lines between z — 2.3 and 
4.1 (Curran ct al. 2011). For comparison, in the Milky 
Way the [Ci]/CO pairs restrict AF/F at the level of 
|AF/F| < 0.4 ppm (Levshakov et al. 2010a) leading to 
a limit on \Aa/a\ < 0.2 ppm since the spatial variations 
in fj, are restricted to |A/i//i| < 0.03 ppm in the disk of 
the Milky Way (Levshakov et al. 2010b; Levshakov et al. 
2010c; Levshakov et al. 2011; Ellingsen et al. 2011). 

The uncertainty of 5 km s" 1 is a statistical error deter- 
mined by the properties of the current observational data, 
such as the S/N ratio, spectral resolution, and spectral line- 
intensity gradients. However, there may be systematic er- 
rors that contribute to the total error budget. In the present 
case, two potential sources of systematic errors are appar- 
ent: (i) velocity shifts caused by observational conditions, 
and (ii) velocity shifts due to non-co-spatial distributions of 
the compared species. The first source is eliminated thanks 
to the unique characteristics of the emission pair [Ci]/CO 
of the z — 5.2 galaxy - narrow linewidths of the subcompo- 
nents and close rest frame-frequencies - which allow us to 
observe them with the same receiver and within the same 
frequency band. In this way, we can avoid systematic un- 
certainties in the velocity scale calibration during long ex- 
posure times. 

We consider possible kinematic segregation between CO 
and C°. In the Milky Way and in nearby galaxies, the C° 
emission is closely associated with that of CO (for refer- 
ences, see Levshakov et al. 2010a; Curran et al. 2011). The 
carbon-bearing species C°, C + , and CO are observed in 
photodissociation regions (PDRs) - neutral regions where 
chemistry and heating are regulated by the far-UV photons 
(Hollcnbach & Tielens 1999). Photons of energy higher than 
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Fig. 3. Examples of the ® 2 (8v) curves for the CO (7-6) and 
[C i] profiles shown in Fig. [5] The curves were generated by 
a bootstrap method as described in Sect. [2j The red line 
shows Q 2 (Sv) of the original data set. 



11.1 eV dissociate CO into atomic carbon and oxygen. Since 
the C° ionization potential of 11.3 eV is quite close to the 
CO dissociation energy, neutral carbon can be quickly ion- 
ized. This suggests that there is a chemical stratification 
of the PDR between the lines C + /C°/CO with increasing 
depth from the surface of the PDR. The observed correla- 
tion between the spatial distributions of C° and CO can be 
explained by clumpy PDR models (e.g., Meixner & Tielens 
1995; Spaans & van Dishoeck 1997; Papadopoulos et al. 
2004). In the observed z = 5.2 galaxy, the derived excitation 
temperature of the higher-J CO lines, T cx = 40K, and the 
gas density nn 2 = 3.5 x 10 3 cm~ 3 (CRR) just corresponds 
to the critical density of n cr « 3 x 10 3 cm~ 3 required to ex- 
cite collisionally the J=2 level of the ground state triplet of 
C° (e.g., Levshakov et al. 2010a). Thus, we can assume that 
in the present case the distributions of CO and C° closely 
trace each other. At z — 5.2, we also observe integrated 
emission over all molecular gas-clouds within the galaxy. 
This naturally leads to an averaging of the random fluc- 
tuations of the line-of-sight velocity components Vco and 
Vjci] and a suppression of the input produced by possible 
deviations from the co-spatial distribution of the compared 
species and the velocity shift between them. 

The uncertainties in the rest-frame frequencies of the 
analyzed transitions can be neglected since they are much 
smaller than the error in the velocity offset, ctav = 5 
km s" 1 (Klein et al. 1998; Muller et al. 2005). 

Thus, we conclude that the limit obtained on the vari- 
ability of AF/F at z — 5.2 is robust and free from signifi- 
cant systematic errors. The error of 5 km s _1 in the veloc- 
ity offset between CO(7-6) and [Ci](2-1) is dominated by 
the noise fluctuations and can be reduced if observations of 
higher S/N are performed. However, even the present con- 
straint on Aa/a = AV/2c — 2±8 ppm casts doubts on the 
dipole model of Webb et al. (2011), which predicts Aa/a = 
-12 ± 3 ppm for the galaxy HLSJ091828.6+514223. 
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Fig. 4. Bootstrap histogram (probability density function, 
p.d.f.) of the velocity offset AV between the CO(7-6) and 
[Ci](2-1) profiles shown in Fig. [2] Gaussian fit corresponds 
to AV = 0.8 ± 4.6 km s" 1 (la C.L.). 
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